1. IntroductionIn conventional quantum mechanics, a Hamiltonian of a quantum system must be represented by a Hermitian operator, which ensures that not only the eigenvalues of the Hamiltonian are real, but also the time evolution of the quantum system is unitary. However, in the last two decades, non-Hermitian (NH) systems have received considerable attention. Many theories of NH systems with real and complex spectra have been investigated. In 1998, Bender et al. proposed that parity–time reversal (
) symmetric NH Hamiltonians which are invariant under combined space and time reversal still have real and positive energy spectra,[1] proved that the time evolution of
-symmetric Hamiltonians is unitary, and redefined an inner product whose associated norm is positive-definite,[2] and demonstrated that the evolutionary time of
-symmetric NH systems can even be made arbitrarily small without violating the time–energy uncertainty principle.[3,4] Assis et al. proved that the phenomenon that the passage time needed for the evolution can be made arbitrarily small in the quantum brachistochrone problem of
-symmetric systems can also be obtained for NH Hamiltonians in which the
-symmetry is completely broken.[5] Günther et al. demonstrated that the
-symmetric quantum brachistochrone problem can be reanalyzed as a quantum system consisting of a
-symmetric NH component and a purely Hermitian component simultaneously,[6] and proposed that the quantum mechanical brachistochrone system with a
-symmetric Hamiltonian can be Naimark-dilated and reinterpreted as a subsystem of a larger conventional quantum mechanics system in a higher-dimensional Hilbert space governed by a Hermitian Hamiltonian.[7] Kawabata et al. proposed the Naimark extension for quantum measurement by adding an ancilla (a measuring apparatus) and extending the Hilbert space, and found that the complete information retrieval from the environment can be achieved in the
-symmetry unbroken phase, whereas no information can be retrieved in the
-symmetry broken phase.[8] Since the quantum brachistochrone problem appeared in NH systems, optimal-speed evolutions of NH systems have been widely investigated.[9–11] Moreover, the framework for the NH formalism of Hamiltonians has been proposed by Brody et al.[12] and Sergi et al.[13] Based on the previous works, many quantum properties and quantum effects in NH systems have been widely studied.[14–20] These research results show that NH Hamiltonians are useful in theoretical work, and they are also regarded as effective mathematical tools for studying quantum properties of open quantum systems in quantum optics.[21–24]
On the other hand, the quantum speed limit time (QSLT) originates from the Heisenberg uncertainty relation for energy and time, it is conventionally known as the minimum evolutionary time between two distinguishable states of a quantum system, and it becomes a key factor in characterizing the maximum evolutionary speed of quantum systems. For closed quantum systems with unitary time evolution, a unified lower bound of the QSLT is obtained by the Mandelstam–Tamm (MT)-type bound τQSL = πħ/(2ΔE)[25] and the Margolus–Levitin (ML)-type bound τQSL = πħ/(2E),[26] where ΔE is the variance of energy of the initial state and E is the mean energy with respect to the ground state. Both the MT-type and the ML-type bounds are attainable in closed quantum systems for initial pure states. According to Refs. [4,27], although
-symmetric NH Hamiltonians with real eigenvalues are not Hermitian in the Dirac sense, they do have entirely real spectra and give rise to unitary time evolution. Hence, the authors in Ref. [28] considered that both the MT-type bound as well as the ML-type bound for time-dependent generators of
-symmetric NH systems remain valid. However, it should be noted that
-symmetric Hamiltonians with real eigenvalues are a special class of NH Hamiltonians. For an NH system which is the spontaneous
-symmetry breaking, the time evolution is no longer unitary.[29] Fortunately, Campo et al. considered NH systems for the quantum speed limit (QSL), and proved that the bound to the speed of evolution under NH Hamiltonians still holds.[30] Deffner et al. expressed the QSLT in terms of the operator norm of the non-unitary generator of the dynamics.[31] However, the unified lower bound introduced in Ref. [31] is applicable for a given driving time for pure initial states, it is not feasible for mixed initial states. Later, Zhang et al. proposed a generic bound on the evolutionary time of quantum systems with non-unitary time evolution by using the relative purity, which is applicable to both mixed and pure initial states.[32] In research field of QSLT, to effectively reduce the evolutionary time as well as accelerate the speed of quantum evolution is one of our strong expectations. Although a lot of works have investigated the quantum Brachistochrone problem in NH systems, this paper is aimed to investigate the evolution time in NH systems from the perspective of the QSLT. With the help of numerical calculations, we demonstrate that the non-Hermiticity in Hamiltonians can notably reduce the QSLT of the system, even decrease to arbitrary small one.
This paper is organized as follows. Firstly, definitions of the QSLT are briefly reviewed. Secondly, physical model and NH dynamics are introduced. Thirdly, QSLT of the NH system is investigated with the help of numerical calculations. Finally, conclusion is given.
2. Definitions of quantum speed limit timeIn this section, we briefly review the definitions of the QSLT. Deffner and Lutz[31] derived a unified lower bound of the QSLT which is determined by an initial state ρ0 = |ψ0⟩〈ψ0| and its target state ρτD, and governed by arbitrary time-dependent non-unitary equation of the form
. With the help of the von Neumann trace inequality and the Cauchy–Schwarz inequality, the QSLT is obtained as
![](cpb_29_3_030304/cpb_29_3_030304_eqn1.gif)
where
τD is the actual driving time,
![](cpb_29_3_030304/cpb_29_3_030304_ieqn16.gif)
, and
![](cpb_29_3_030304/cpb_29_3_030304_ieqn17.gif)
denotes the Schatten
p norm,
σ1, …,
σn are the singular values of
A. ||
A||
1 =
∑i σi is the trace norm,
![](cpb_29_3_030304/cpb_29_3_030304_ieqn18.gif)
is the Hilbert–Schmidt norm, and ||
A||
∞ =
σmax is the operator norm which is given by the largest singular value. Because of the relationship ||
A||
∞ ≤ ||
A||
2 ≤ ||
A||
1, the ML-type bound based on the operator norm (
p = ∞) provides the sharpest bound on the QSLT. And
![](cpb_29_3_030304/cpb_29_3_030304_ieqn19.gif)
is the Bures angle between the initial state
ρ0 and the target state
ρτD. However, equation (
1) is not feasible for mixed initial states. Fortunately, Zhang
et al. in Ref. [
32] proposed a unified lower bound of the QSLT for an arbitrary mixed state
ρτ to its target state
ρτ + τD based on the relative purity:
![](cpb_29_3_030304/cpb_29_3_030304_eqn2.gif)
where
![](cpb_29_3_030304/cpb_29_3_030304_ieqn20.gif)
,
σi and
ρi are the singular values of
Ltρt and the mixed initial state
ρτ, respectively. And
![](cpb_29_3_030304/cpb_29_3_030304_ieqn21.gif)
denotes the relative purity between the initial state
ρτ and the final state
ρτ + τD with the driving time
τD. For a pure initial state
ρτ = 0 = |
ψ0⟩〈
ψ0|, the singular value
ρi =
δi,1, then
![](cpb_29_3_030304/cpb_29_3_030304_ieqn22.gif)
, and equation (
2) can be simplified as
![](cpb_29_3_030304/cpb_29_3_030304_eqn3.gif)
Equation (
3) indicates that, with regard to a pure initial state, the expression given by Eq. (
2) can recover to the unified lower bound of the QSLT obtained by Eq. (
1). When
τQSL =
τD, it means that the evolution is already along the fastest path and does not possess potential capacity for further quantum speed-up. However,
τQSL <
τD indicates that the further acceleration might occur, and the shorter
τQSL, the greater the capacity for potential speed-up will be. And 1/
τQSL defines a natural notion of the speed of quantum evolution. Because two bounds of the QSLT can be saturated (
i.e., the actual optimal evolution time could equal the speed limit time), reducing
τQSL would lead to an acceleration of the quantum evolution. Especially,
τQSL = 0 can be interpreted as two different situations: for two identical states (
i.e., the initial state is the same as its target state),
τQSL = 0 indicates that the quantum evolutionary speed tends to zero-speed, but for two distinguishable quantum states (
i.e., the initial state is different to its target state),
τQSL = 0 represents that the quantum evolutionary speed becomes infinite speed.
3. Physical model and non-Hermitian dynamicsWe all know that the NH approach is regarded as one of available methods to describe properties of open quantum systems. Sergi et al.[13] assumed that in absence of any interaction with the environment, the two-level system is free to make transitions between its two energy levels. Such a situation is modeled by the Hermitian Hamiltonian H+. In order to formulate the open system dynamics of the model, they introduced a general anti-Hermitian Hamiltonian H−. Namely, NH Hamiltonians (
) always can be decomposed into Hermitian and anti-Hermitian parts as HnH = H+ + H− with
and H− = −iΓ, where Γ = Γ† is usually regarded as the decay rate operator. We choose a special scheme which has been realized in both classical experiments[33–35] and quantum experiments[36] to describe the NH model for which gain and loss of energy or amplitude are present (ħ = 1):
![](cpb_29_3_030304/cpb_29_3_030304_eqn4.gif)
where
ω and
γ are assumed to be real-valued,
γ represents the NH parameter, and
σα (
α =
x,
z) are Pauli matrixes. And the Hamiltonian of the NH model is given as
![](cpb_29_3_030304/cpb_29_3_030304_eqn5.gif)
where we denoted
Δ =
γ/
ω, and
ω−1 is an energy scale, namely
Δ has the same sign as the NH parameter
γ. And the eigenvalues of
HnH are
![](cpb_29_3_030304/cpb_29_3_030304_ieqn25.gif)
. It is easy to determine that when
γ ∈ (−1,0) ∪ (0,1),
HnH in the
![](cpb_29_3_030304/cpb_29_3_030304_ieqn26.gif)
-symmetric phase with real eigenvalues, while
γ ∈ (−∞,−1) ∪ (1,∞),
HnH in the
![](cpb_29_3_030304/cpb_29_3_030304_ieqn27.gif)
-symmetry broken phase with complex eigenvalues. And
γ = ±1 are usually considered as exceptional points (EPs) where eigenvalues switch from real values to complex values.
[37] It should be noted that
![](cpb_29_3_030304/cpb_29_3_030304_ieqn28.gif)
-symmetry breaking is usually associated with the EP; in contrast, in the system with infinite site, the
![](cpb_29_3_030304/cpb_29_3_030304_ieqn29.gif)
-symmetry breaking is associated with spectral singularity rather than the EP.
[38,39] When
γ = 0, the NH model degrades into a coherent Rabi oscillation with coupling parameter
ω = 1. The evolution equation of NH systems with an initially pure state is a complex extension of Schrödinger equation. For a mixed state, the evolution equation of NH systems refers to the covariance equation which is a complex extension of the von Neumann equation
[12]
![](cpb_29_3_030304/cpb_29_3_030304_eqn6.gif)
where [ , ] and { , } represent the commutator and the anti-commutator, respectively. In general, due to the non-Hermiticity of
HnH, the NH dynamics is non-unitary, and
Ωt in Eq. (
6) is a non-normalized density operator. Therefore, for making sure that the density matrix is trace-preserving, the following renormalization process is required:
![](cpb_29_3_030304/cpb_29_3_030304_eqn7.gif)
And then the norm-preserving evolution equation generated by an NH Hamiltonian for the normalized density operator
ρt is given as
![](cpb_29_3_030304/cpb_29_3_030304_eqn8.gif)
It should be noted that the solution of this evolution equation in Eq. (
8) also can be expressed in the form
![](cpb_29_3_030304/cpb_29_3_030304_eqn9.gif)
where
ρ0 is the initial state,
UnH = exp(−i
HnHt) is still a non-unitary time evolution operator. In conventional quantum mechanics, the NH dynamics is intrinsically non-unitary both in the
![](cpb_29_3_030304/cpb_29_3_030304_ieqn30.gif)
-symmetry unbroken and broken phases.
[33] However, the Hermitian dynamics can appear in non-Hermitian systems.
[40]By solving the evolution equation Eq. (8), matrix elements
(k,l = 1,2) of the normalized final state ρt are given as
![](cpb_29_3_030304/cpb_29_3_030304_eqn10.gif)
where
![](cpb_29_3_030304/cpb_29_3_030304_ieqn32.gif)
(
i,
j = 1,2) are elements of the initial state
ρ0, we have denoted
![](cpb_29_3_030304/cpb_29_3_030304_ieqn33.gif)
and
![](cpb_29_3_030304/cpb_29_3_030304_ieqn34.gif)
.
4. Quantum speed limit time of non-Hermitian modelAccording to definitions of the QSLT given by Eqs. (1) and (2), we consider two cases of different initial states ρ0.
4.1. Pure initial state caseWe firstly examine the QSLT of the NH two-level system with a pure initial state |ψ0⟩ = |1⟩, and equation (1) can be simplified as
![](cpb_29_3_030304/cpb_29_3_030304_eqn11.gif)
where
![](cpb_29_3_030304/cpb_29_3_030304_ieqn35.gif)
represents the population of the excited state |1⟩ at time
τD, and
σmax is the largest singular value of
Ltρt.
In Fig. 1, we depict the QSLT represented by τQSL (the black-solid curve) as a function of the NH parameter represented by Δ, and the actual driving time τD = 1. According to Eq. (5) and its explanation, when Δ ∈ (−1,1), HnH is a
-symmetric NH Hamiltonian, while Δ∈(−∞,−1)∪(1,∞), HnH is a
-symmetry broken NH Hamiltonian, and Δ = ± 1 are usually considered as EPs. On the basis of definitions of the QSLT, and the actual optimal evolution time could equal the speed limit time, we know that increasing τQSL would lead to a deceleration of the quantum evolution, while decreasing τQSL would result in an acceleration, and τQSL = 0 can be interpreted as zero-speed or infinite speed. Besides, from Eq. (11), one can see that τQSL is related to p(τD), and p(τD) can reflect evolutionary efficiency of the quantum state. In order to explore the internal mechanism of the speed of the quantum evolution, we also plot p(τD) (the red-dashed curve) as a function of Δ. In Fig. 1, we adjust Δ from negative values to positive values. For negative and small region Δ ∈ (−15,−3), p(τD) slightly decreases, and τQSL increases correspondingly. Especially, when Δ ∈ (−3,1.5), p(τD) quickly decreases, and τQSL firstly increases, and then decreases. For positive and large region Δ∈(1.5,15), p(τD) slightly reverts, and finally decreases to zero, τQSL firstly decreases and then increases to a constant. That is to say, when Δ is in negative and small region, although the value of τQSL is smaller, p(τD) barely changes, which means that the evolution of the quantum state is faster but low efficient. In the region which p(τD) rapidly and efficiently changes, the quantum evolution experiences a process from deceleration to acceleration, and then τQSL gradually tends to a constant speed with Δ increases continuously. It is obvious that, with regard to the NH parameter, the excited state population has a negative relation with the QSLT. The speed-up of NH systems is mostly because that the major difference between conventional and NH quantum mechanics is the definition of the inner product. Because the Hilbert-space metric depends on the Hamiltonian, the geometry of Hilbert space of the NH quantum theory has to be modified. Hence, a pair of states is orthogonal under the standard inner product of the Hermitian quantum theory, but is no longer orthogonal in the NH quantum theory. As a consequence, the counterpart of Bures angle in Eq. (1) under the NH quantum theory must be changed, which is possible that an alternative complex pathway from a state to its target state can be made shorter.
We secondly explore effects of the non-Hermiticity on the QSLT of the two-level system in the whole dynamical process by using the definition given by Eq. (2). We also start from the excited state |1⟩ and discuss the QSLT τQSL from an arbitrary state ρτ to its target state ρτ + τD, and
represents the population of the excited state |1⟩ at the time τ. In Fig. 2, we plot the QSLT τQSL (black-solid curves) and the excited state population p(τ) (red-dashed curves) as functions of the initial time parameter τ for different NH parameter Δ, and the actual driving time τD = 1. According to Eq. (5) and its explanation, we know that HnH with the NH parameter Δ = 0.4 or 0.9 in Fig. 2(a) or Fig. 2(b) satisfies
-symmetric structure. In Figs. 2(a) and 2(b), τQSL periodically decreases and increases, which implies that the quantum evolution exists periodical speed-up and speed-down. We also find that with Δ increases from 0.4 to 0.9, the effect of acceleration and deceleration becomes more obvious. It is means that, in the
-symmetric phase, the larger the NH parameter Δ is, the more obvious the effect of speed-down and speed-up is, the smaller value the QSLT can be achieved, which corresponds to a faster speed of the quantum evolution. And p(τ) periodically evolves and has the same period with τQSL. While Δ = 1.1 or 2.5 in the
-symmetry broken phase in Fig. 2(c) or Fig. 2(d), the quantum evolution is aperiodic. From Figs. 2(c) and 2(d), we can see that τQSL slightly increases at the beginning, and gradually decreases to zero, which implies that the quantum evolution exists a speed-down and a speed-up. And p(τ) firstly decreases and then increases, finally reaches at a stable value. The zero QSLT combined with a stable excited state population indicates that the state of the NH system finally evolves into a steady state and the quantum evolutionary speed finally tends to zero-speed. In the
-symmetry broken phase, the larger the NH parameter Δ is, the faster the evolution of the system to the steady state is, the smaller the stable value of the excited state population p(τ) is, which corresponds to faster and more effective quantum evolution.
4.2. Mixed initial state caseIn this subsection, we examine the QSLT of the NH two-level system with a mixed initial state
![](cpb_29_3_030304/cpb_29_3_030304_eqn12.gif)
where
p is a constant parameter. For clearly demonstrating the quantum evolution process, we also consider the trace distance
[41] which is a measure of the distinguishability between two quantum states
ρ1 and
ρ2
![](cpb_29_3_030304/cpb_29_3_030304_eqn13.gif)
with
![](cpb_29_3_030304/cpb_29_3_030304_ieqn43.gif)
. The trace distance of two distinguishable states satisfies the inequality 0 <
D(
ρ1,
ρ2) ≤ 1,
D(
ρ1,
ρ2) = 0 for two identical states
ρ1 =
ρ2, and
D(
ρ1,
ρ2) = 1 for two orthogonal states
ρ1 ρ2 = 0.
Effects of the non-Hermiticity on the QSLT of the two-level system in the whole dynamical process can be studied by using Eq. (2) when the initial state takes the form of a general mixed state given by Eq. (12). In Fig. 3, we depict the QSLT τQSL (black-solid curves) and the trace distance D(ρτ,ρτ + τD) (blue-dotted curves) as functions of the initial time parameter τ for different NH parameter Δ. According to Eq. (5) and its explanation, we know that HnH with Δ = 0.6 or 0.9 in Fig. 3(a) or Fig. 3(b) satisfies
-symmetric structure, and the quantum evolution is periodical. We find that, in Fig. 3(a), Δ = 0.6, the quantum evolution exhibits straightforward acceleration and deceleration in one period. While in Fig. 3(b), Δ = 0.9, τQSL tends to zero twice in one period, and the quantum evolution experiences repeated speed-down and speed-up. It is worth pointing out that when τQSL reduces to a minimum, D(ρτ,ρτ + τD) is nonzero. The zero QSLT combined with the nonzero trace distance indicates that the quantum state of the NH system still evolves, but evolutionary time is zero, which means that the quantum evolutionary speed tends to infinite speed. That is to say, the evolutionary time of NH systems do not have a nonzero lower bound, which is a remarkable difference comparing with the conventional quantum theory. In order to give NH quantum systems a meaning in the conventional quantum mechanics, the authors in Ref. [8] embedded NH systems into a larger Hermitian system, and proposed the Naimark extension for quantum measurement by adding an ancilla (a measuring apparatus) and extending the Hilbert space, any non-unitary dynamics can be understood as a unitary dynamics of the entire system followed by quantum measurement acting on the ancilla. Moreover, the authors in Ref. [7] proposed that the embedding of NH system into a higher-dimensional Hilbert space can be deemed as a strengthening of the wormhole analogy introduced in Ref. [3], and the authors in Ref. [3] noted that it is possible to create a wormhole-like effect in the Hilbert space to explain why the transformation between a pair of orthogonal states (under the standard inner product in Hermitian quantum theory) can be made in arbitrarily small time. This is because that for a quantum system described by NH Hamiltonians the alternative complex pathway from a state to its orthogonal state can be made arbitrarily short. The mechanism described here is similar to that in general relativity in which the alternative distance between two widely separated space–time points can be made small if they are connected by a wormhole. Besides, an exceptional point (EP) of the NH Hamiltonian has only one (geometric) eigenvector, since both the eigenvalues and the corresponding eigenstates of the NH Hamiltonian coalesce at EPs, in Ref. [11] the authors found that some NH Hamiltonians correspond to an NH degeneracy called EPs, and showed that any state evolution can be generated solely by such NH degeneracies yielding an EPs-driven evolution which minimizes the Hilbert–Schmidt norm of the matrix of NH Hamiltonians. And in Ref. [8] the authors proposed that the EPs plays a role of the critical points around which many physical quantities such as the recurrence time and the distinguishability show power-law behavior. These findings may find novel applications to quantum control. Hence, in the following, we choose some typical and representative cases (i.e., EPs) to expound our work. In Figs. 3(c) and 3(d), Δ = 1 and −1 are EPs of
-symmetric structure in the NH model, and the quantum evolution is aperiodic. We can also find that when τQSL decreases to zero, D(ρτ,ρτ + τD) is still nonzero, which is the same as the phenomenon that appears in Figs. 3(a) and 3(b). The nonzero trace distance and the zero QSLT indicates that the quantum evolutionary speed tends to infinite speed. In addition, we also consider others situations of Δ. Our numerical calculations also show that the quantum evolutionary speed of distinguishable states of the NH two-level system can tend to infinite speed. It is well known that the surprising result in Refs. [3,5] showed that the evolutionary time of NH systems can be made arbitrary small without violating the time–energy uncertainty principle, and we demonstrated the result by using QSLT in this work.
5. ConclusionIn this paper, we considered an NH two-level system in the presence of gain and loss of energy or amplitude, and studied the QSLT of the NH system in regard to two cases of pure and mixed initial states. In the pure initial state, the quantum evolution of the system can be effectively accelerated by adjusting the NH parameter. While in the mixed initial state, with respect to two distinguishable states of the NH system, the evolutionary time does not have a nonzero lower bound, as well as the QSLT can be arbitrarily small even be zero. A probable physical explanation for why the NH systems allow for faster evolutions is that the major difference between conventional and NH quantum mechanics is the definition of the inner product. Hence, the geometry of Hilbert space of the NH quantum theory has to be modified. It is easy to note that, the modification of NH Hilbert space has to be implemented, irrespective of whether two distinguishable states of NH systems are orthogonal or not. That is to say, under the NH quantum theory, by adjusting the non-Hermiticity of the system, the alternative complex pathway of such two distinguishable states can be made in arbitrarily short, and the evolutionary time of the states also can be arbitrarily small without violating the time–energy uncertainty principle.